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Abstract 


We examine whether the inverted hierarchical model of neutrinos is compat¬ 
ible with the explanation of the large mixing angle (LMA)MSW solution of 
the solar neutrino problem. The left-handed Majorana neutrino mass matrix 
for the inverted hierarchical model, is generated through the seesaw mech¬ 
anism using the diagonal form of the Dirac neutrino mass matrix and the 
non-diagonal texture of the right-handed Majorana mass matrix. In a model 
independent way we construct a specific form of the charged lepton mass 
matrix having a special structure in 1-2 block, which contribution to the 
leptonic mixing (MNS) matrix leads to the predictions sin 1 2 29\2 = 0.8517, 
sin 2 20 2 3 = 0.9494 and 14,3 = 0.159 at the unification scale. These predic¬ 
tions are found to be consistent with the LMA MSW solution of the solar 
neutrino problem. The inverted hierarchical model is also found to be stable 
against the quantum radiative corrections in the MSSM. A numerical anal¬ 
ysis of the renormalisation group equations (RGEs) in the MSSM shows a 
mild decrease of the mixing angles with the decrease of energy scale and the 
corresponding values of the neutrino mixings at the top-quark mass scale are 
sin 2 20i2 = 0.8472, sin 2 29 2 3 = 0.9399 and |14 3 | = 0.1509 respectively. 
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1 Introduction 


Neutrino physics is one of the fast developing areas of particle physics. 
The recent Super-Kamiokande experimental results on both solar[1] and at¬ 
mospheric [2] neutrino oscillations support the approximate bimaximal mix¬ 
ings. Though these results favour the large mixing angle (LMA) MSW so¬ 
lution with active neutrinos, such interpretation is not beyond doubt at this 
stage[3,4,5]. We also have little idea about the pattern of the neutrino mass 
spectrum whether it is hierachical or inverted hierarchical, and both possi¬ 
bilities are consistent with the neutrino oscillation explanations of the at¬ 
mospheric and solar neutrino deficits[5,6]. The data from the long baseline 
experiment using a Neutrino factory will be able to confirm the actual pattern 
of the neutrino masses in the near future [7]. 

In the theoretical front the hierarchical model of neutrino masses and 
its generation have been widely studied and found to be consistent with 
the explanation of the LMA MSW solar neutrino solution[8,9]. However, 
the inverted hierarchical model of neutrino masses generally predicts the 
maximal mixing angles 9 and 9 23 close to 45°, and are suitable for the 
explanation of the vacuum oscillation (VO) solution of the solar neutrino 
oscillation[6,10] and the atmospheric neutrino oscillation data. The atmo¬ 
spheric data gives the lower bound at sin 2 29 2 3 > 0.88 and the best-fit value 
at Am 23 = 3 x lCT 2 eV 2 . It is quite obvious that the prediction from the 
inverted hierarchical model fails to explain the LMA MSW solution which 
has upper experimental limit[4,6] sin 2 29\2 < 0.988 at 95% C.L., and the 
best-fit values sin 2 29\2 = 0.8163 and Am 2 2 = 4.2 x 10 -5 eV 2 . Combining 
LMA MSW solution and atmospheric data, the best-fit value of the mass 
splitting parameter is obtained[6] as £ = Am 2 2 /Am 23 = 0.014. It has been 
argued [6,10] that the diagonalisation of the charged lepton mass matrix 
cannot give a significant contribution to 9 u Vi needed for the explanation of 
the LMA MSW solution. On such ground the inverted hierarchical models 
are assumed to be inconsistent with the LMA MSW solution. An attempt 
was made to explain the LMA MSW solution from the inverted hierarchical 
model by considering two types of charged lepton mass matrices[11] and was 
partially successful. We are interested here to make further investigations in 
this paper whether the inverted hierarchical model gives an acceptable LMA 
MSW solution when we include the contribution from the diagonalisation of 
the charged lepton mass matrix having special form in the 1-2 block, to the 
leptonic mixing matrix. 
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The paper is organised as follows. In section 2 we outline the seesaw 
mechanism for the generation of the neutrino mass matrix which leads to the 
inverted hierarchical mass pattern, and the contraction of the charged lepton 
mass matrix suitable for the LMA MSW solution. In section 3 we describe 
briefly the procedure for the analysis of the renormalisation group equations 
(RGEs) within the minimal supersymmetric standard model (MSSM). This 
is followed by a summary and discussion in section 4. 


2 Generation of the inverted hierarchical neu¬ 
trino mass matrix and the charged lepton 
mass matrix 

The inverted hierarchical model of neutrinos has its origin from the low 
energy non-seesaw models [12], e.g., the Zee-type of model using a singly 
charged singlet scalar held and also the models with an approximate con¬ 
served L e — — L t lepton number. However, it is also possible to generate 

the inverted hierarchical model through the seesaw mechanism at high energy 
scale within the framework of the grand unified theories with a chiral U(l) 
family symmetry[10,ll]. In a model independent way, we consider the Dirac 
neutrino mass matrix uilr and the non-diagonal form of the right-handed 
Majorana mass matrix M R in the seesaw formula [13] given by 

m LL = m LR M R x m T LR (1) 

where rrijj J is the left-handed Majorana mass matrix. The leptonic mixing 
matrix known as the MNS mixing matrix [14] is defined by 

Vmns = KlKl (2) 

where 1and V,l are obtained from the diagonalisation of the charged 
lepton and uill as 

my = V ei m e Vj R 

'"'ll 1 = V„ L m LL Vj L (3) 

If the charged lepton mass matrix is diagonal, the MNS matrix (2) is simply 
given by 

Vmns = Vj L (4) 
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We can always express will in the basis where the charged lepton mass matrix 
is diagonal, 

™'ll = Ve L m LL V} L , 

Idiag _ y, , yl\ 
m LL ~ V uL m LL V eLi 

Vmns = vX (5) 

From the above expressions we can calculate the experimentally deter¬ 
mined quntities as follows: 

(i) the neutrino mass splitting parameter, £ = ^"'l 2 

(ii) the atmospheric mixing angle, S at = sin 2 2d 23 = 4| V^, 3 1 2 (1 — | V^, 3 1 2 ) 

(iii) the solar mixing angle, S so i = sin 2 2 #i 2 = 41 V^ 2 1 2 1 | 2 

(iv) the CHOOZ angle, S c = 4|F4 3 | 2 (1 - |F4 3 | 2 ) or simply \V eZ \. 

The Vfi where / = r, /i, e and * = 1,2, 3, are the elements of the MNS mixing 
matrix. 

First, we consider the diagonal form of the charged lepton mass matrix 
m e given by 

/A 6 0 0 \ 

rn e — 0 A 2 0 m T (6) 

V o 0 1 ) 

where the Wolfenstein parameter[15] A = 0.22 and the ratios of the charged 

lepton masses are m T : : m e = 1 : A 2 : A 6 respectively. From Eq.(6) we get 

V e L = 1 and Vmns = VvL as i n Eq. (4). Again we consider the diagonal form 
of the Dirac neutrino mass matrix wilr as the up-quark mass matrix[16], 

/A 8 0 0 A 

m L R = 0 A 4 0 I m t (7) 

v o o i; 

where the up-quark masses are in the ratios[17] m t : m c : m u = 1 : A 4 : A 8 . 
Now, the proper choice of the elements in M R , enables us to generate the 
inverted hierarchical neutrino mass matrix. We present here the following 
examples [18]: 


Example (a) 


( 0 1 1 

will = 1 A 3 0 

V 1 0 A 3 
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with the choice 


Eq.(8) yields 


/ -A 


22 


Mr = 


V A 


A 15 

11 


A 15 A 11 

A 8 -A 4 
-A 4 1 


vr 


Vmns 


K f L = 


/ 0.70577 
-0.50094 
V -0.50094 


0.70844 -1.11 x 10” 16 
0.49906 -0.70711 

0.49906 0.70711 


vr 


(9) 


and the neutrino mass eigenvalues m, = (1.4195,1.4089,0.0105)mo, i = 
1,2,3. This gives the mass splitting parameter £ = Araf 2 /Am| 3 = 0.014, 
and the mixing angles sin 2 20 12 = 0.9999, sin 2 20 2 a ~ 1.00, V e3 = 0.0. 


Example (b) 


m LL 


( 0 1 
1 —(A 3 — A 4 ) / 2 
\1 —(A 3 + A 4 )/2 


—(A 3 + A 4 )/2 mo, 
-(A 3 - A 4 )/2 j 


( 10 ) 


with the choice 


/ 

A 23 

A 16 

A 12 


A 16 

A 8 

-A 4 

V 

A 12 

-A 4 

1 


leading to rry = (1.4195,1.4089, 0.00239)m o , ^ 
sin 2 2023 ~ 1-00, |Ks| = 0.0. 


vr 

0.0151, sin 2 20 12 


0.9999, 


Example (c) 


m LL 


I A 3 1 1 \ 

1 A 4 /2 —A 4 /2 m 0 , 

V 1 —A 4 /2 A 4 /2 j 


with the choice 


Mr = 


( 0 

A 16 


A 16 

A 8 

— (A 4 + A 12 ) 


A 12 


— (A 4 + A 12 ) 
1 


vr 


( 11 ) 
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leading to m t = (1.4195,1.4089, 0.002343)m o , $, = 0.015, sin 2 20 12 = 0.9999, 

sin 2 2023 ~ 1-00, \Va\ = 0.0. 

In the above results the Vmns obtained from the uill alone fails to explain 
the LMA MSW solution, and any small deviation in itill will hardly affect 
sin 2 2di2 [6,10]. The last hope is that there could be a significant contribution 
to 0 i 2 from the diagonalisation of the charged lepton mass matrix having 
special structure in 1-2 block [11], We wish to examine how 9 so i = (0 32 — 0 32 ) 
can resolve the LMA MSW solar neutrino mixing scenario [11]. 

We parametrise the charged leptonic mixing V e L by the following three 
rotations[19,20] 

VeL = A23A13A12 


/1 

0 

0 ^ 

( C13 

0 

S13 A 

( 

C12 

S12 

0 A 

0 

C23 

S23 1 

0 

1 

0 


-S12 

C12 

0 ( 12 ) 

Vo 

—S23 

C23 / 

V -«13 

0 

Cl 3 j 

V 

0 

0 

1 / 


where 

to 


sin 0^- and c l3 = cos 0?-. Putting 0 33 = 0| 3 


V eL = 


( C 12 


-S12 

V 0 


S 12 

C 12 

0 


0 

0 

1 


0, Eq.(12) reduces 


(13) 


This gives a special form in the 1-2 block. We then reconstruct [19] the 
charged lepton mass matrix using Eq.(13) from the relation, 


m e = V} L mf aa V eR 


f X 6 cf 2 + A 2 S 2 2 A 6 Ci 2 Sl2 - A 2 Ci 2 Sl2 0 A 

= A 6 Ci 2 Sl2 - A 2 Ci 2 Sl2 A 6 c 2 2 + A 2 s 2 2 0 (14) 

Vo 01 ) 

For a specific choice of d\ 2 = 13°, and A = 0.22, Eq.(14) leads to 


rn e 


( 0.00256 -0.01058 0 A 

-0.01058 0.04596 0 

V 0 Oly/ 


(15) 


which has a special form in the 1-2 block. The diagonalisation of m e in 
Eq.(15) gives 


V eL = 


-0.97439 

-0.22488 

0 


-0.22488 0 
0.97439 0 

0 1 


(16) 
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which is now completely unitary. The corresponding eigenvalues of the 
charged lepton mass matrix are given by 

mf a9 = (1.182 x 10' 4 ,4.84 x 10" 2 ,1.0)m T (17) 

which gives almost correct physical mass ratios[17] m e : m M : m T = A 6 : A 2 : 1. 
The MNS mixing matrix (2) is now calculated, using Eqs.(9) and (16), as 


Vmns 


v eL v} 


vL 


I -0.5750 
-0.6468 
V -0.50094 


-0.8025 

0.32696 

0.49906 


0.1590 \ 
0.6890 
0.70711 ) 


(18) 


This leads to the mixing angles sin 2 29\ 2 = 0.8517, sm 2 20 2 3 = 0.9494, and 
V e3 = 0.159, and these predictions are consistent with the explanation of 
LMA MSW solution. The possible choice of 9f 2 = 14° in Eq.(14) also leads 
to the predictions of sm 2 2# 12 = 0.8298, sin 2 29 23 = 0.9415, and \V e3 \ = 
0.1710 while maintaining the good prediction of the ratios of the charged 
lepton masses. However its \V e3 \ is above the CHOOZ and PALO VERDE 
experimental constraint [21] of \V e3 \ < 0.16. 

Taking the first result with 9\ 2 = 13°, the left-handed neutrino mass m' LL 
in the basis where the charged lepton mass matrix is diagonal(5), is now 
expressed for our convenience, as 


/ 0.437972 
-0.897698 
V -0.973193 


-0.897698 -0.973193 
-0.443296 -0.230068 
-0.230068 -0.005324 


m 0 


(19) 


where Vmns — V„ L is same as earlier given in Eq.(18), and the neutrino mass 
eigenvalues are 

rn t = (1.4196,1.4089,4.234 x 10~ 7 ) m 0 ; i = 1,2,3 
which give the mass splitting parameter, £ = Am 2 2 /Am 23 = 0.014. 


3 Renormalisation effects in MSSM 

It is desirable to inspect how the values of sin 2 29i2, sin 2 29 2 3, \V e3 \ and ^ 
evaluated at the unification scale where the seesaw mechanism is operative, 
respond to the renormalisation group analysis on running from higher scale 
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(M u = 2 x 10 16 GeU) down to the top quark mass scale (// = m t ) [22], We 
consider the renormalisation group equations (RGEs) for the three gauge 
couplings (<?i, <?2, < 73 ) and the third family Yukawa couplings ( h t ,hb,h T ) in 
the minimal supersymmetric standard model (MSSM) in the standard fash- 
ion[23]. At high scale /a = 2 x 10 16 GeU, we assume the unification of gauge 
couplings as well as third generation Yukawa couplings for large tan/3 [23]. 
We choose the input a 2 = (5/3)aq = a 3 = 1/24, and h t = hi , = h T = 0.7 
corresponding to large tan/3 = v u /vd- 

We express nIII in terms of K, the coefficient of the dimension five neu¬ 
trino mass operator[24,25,26,27] in a scale-dependent manner, 

m LL (t) =v 2 u (t)K(t) (20) 

where t = ln(/a) and v u (t) is the scale-dependent[27] vacuum expectation 
value (VEV) v u = vosin/3, Vo = 174 GeV. In the basis where the charged 
lepton mass matrix is diagonal, we can write Eq.(20) as [25,27] 

m' LL (t) = vl(t)K\t) (21) 


where K'(t ) is the coefficient of the dimension five neutrino mass operators 
in the basis where the charged lepton mass matrix is diagonal. The evolution 
equations are given by[27] 


dt lnv n(t) 167r2 L20 


3 2 3 9 „, 2 

- 3 h t 


* lnK '® = -ih 


\gl + 6^2 - 6 h 2 t - h 2 r S i: 

5 


hlS, 


The evolution equation of m' LL (t) in Eq.(21) simplihes[27] to 


j t lnm' LL (t) = W 


~Iq9i - - 2 hl + h 2 T Si3 + h%j 


Upon integration from high scale (t u = In M u ) to lower scale (to 
where t 0 < t <t u and t = In /i, we get [27] 


x 


m' iL (t 0 ) = e-To In(to) e^ Ig ^ to) 

(Q rn’ LL12 (t u ) m' LL i 3 (t u )e~ lT{to) 

™>L L2 i(tu) rn' LL22 (t u ) rn' LL23 (t u )e~ IAto) 

V m LL3i(tu)e- IAto) m' LL32 (t u )e- lT{to) vn! LL23 (t u )e~ 2I ^ 


( 


( 22 ) 

( 23 ) 

(24) 
In m t ) 


m 


LLll \ v u 


(25) 
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where 


4(*o) = ~^~2 J 9i(t)dt, * = 1,2, 3; 

1 C^ u o 

If (to) = J t h)(t)dt , f = t,b, T. (26) 

Using the numerical values of I gi (t ) and //(£) at different energy scales t, t 0 < 
t < t u the left-handed Majorana mass matrix m' LL (t ) in Eq.(25) is estimated 
at different energy scales from the value of rn' LL (t u ) given in Eq.(19). At 
each scale the leptonic mixing matrix Vmtvs(^) = V^(i) is calculated and 
this in turn, gives mixing angles sm 2 2d 12 , sm 2 2# 23 and V e3 \. For example, 
at the top-quark mass scale t 0 = In rn t = 5.349, we have calculated I T (t 0 ) = 
0.100317 and the leptonic mixing matrix 


Vmns 


f -0.56962 
-0.67061 
V -0.4752 


-0.80795 

0.35075 

0.47349 


0.15085 \ 
-0.65364 
0.74161 ) 


(27) 


which leads to the low-energy predictions: sin 2 29i2 = 0.8472 and sm 2 2d 23 = 
0.9399. There is a mild reduction from the values estimated at the high 
energy scale t u . This feature shows the compatibility of the inverted hier¬ 
archical model with the explanation of the LMA MSW solution. The pa¬ 
rameter |U e3 | = 0.15085 meets the CHOOZ constraint \V e3 \ < 0.16 [21]. 
The neutrino mass eigevalues at low-energy scale are obtained as m, = 
(1.3533,1.3436,3.8376 x 10 _7 )m o . However the mass splitting parameter 
^ = Am 2 2 /Am 23 = 0.01449 remains almost constant. The running of the 
mixing angles S so i = siri 2 29i2 and S a t = sin 2 2d 23 is shown in Fig.l by the 
solid-line and dotted-line respectively. Both parameters decrease with de¬ 
crease in energy scale t. This is a desirable feature for maintaining the 
stability condition of the inverted hierarchical model. Note that there is 
no exponential term which depends on the top-quark Yukawa coupling in¬ 
tegration in Eq.(25). This differs from the expressions calculated in earlier 
papers[ll,25]. The absence of such term increases the stability criteria of 
the mass matrix at low tan/3 region and this is made possible only when 
we consider the scale-dependent vacuum expectation value in Eq.(22). As 
discussed before, if the CHOOZ constraint is relaxed to \V e3 \ < .2, then we 
would be able to get even lower value of sin 2 29\2 suitable for the explanation 
of the best-fit value of the LMA MSW solution. 




5 10 15 20 25 30 35 40 

t = ln(fi) 

Fig.l. Variations of S so i = sin 2 26\2 and S at = sm 2 20 2 3 with energy scale 
t = ln(/i), which are represented by solid-line and dotted-line, respectively. 

We now briefly discuss the analytic solution for the evolution of the neu¬ 
trino mixings in MSSM. The equation for the evolution of the neutrino mix¬ 
ing angle sin 2 29 can be approximated as two flavour mixing in terms of the 
neutrino mass eigenvalues[26] mi,m 2 , m 3 by a generalised evolution equation, 

167T 2 —— sin 2 2 Oij = —2 sin 2 26L(1 — sin 2 2%)(h 2 — h?) — + — (28) 

dt 3 rrij — rrii 

where i < j and i,j = 1, 2, 3; or e, //, r. For solar mixing angle 9 V2l we have 

167r 2 -^-sin 2 2012 = —2 sin 2 20 12 (1 — sin 2 20 12 )(h 2 — h 2 ) m2 — (29) 

dt A m 2 — mi 
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and for atmospheric mixing angle 6*23 we have 

167 r 2 -^- sin 2 26*23 = —2 sin 2 2023(1 — sin 2 29 23 )^ — h 2 ) — — — (30) 

ere ^ m 3 — 771 2 

Since /r 2 > h 2 > /r 2 and mi > m 2 > m 3 , both 5 so ; and 5 at in Eqs.(29) and 
(30) increase with increasing energy scale. These features are plainly visible 
in the Fig.l. 

A few more comments on our choice of the texture of the charged lepton 
mass matrix are also presented. Here we examine the forms of the texture 
of the charged lepton mass matrix m e and its diagonalisation matrix V e l 
in Eqs.(15)and (16). For simplicity of our analysis we express them in the 
following approximate forms 


and 


rn e ~ 


( A 4 
-A 3 

V 0 


-A 3 0 A 
A 2 0 
0 1 ) 


V eL ~ 


-1 -A 0 A 
-A 1 0 

0 0 1 / 


(31) 


(32) 


It is interesting to note that the position of the zeros in the mass matrix in 
Eq.(31) have the same structure with those of lepton mass matrix obtained 
by Ibanez and Ross[29] in the gauge theory of the standard model with an 
horizontal U( 1) gauge factor. Such form of the texture of the charged lepton 
mass matrix is also proposed by Georgi and Jarlskog[30] in 577(5) model, 
and this can be realised in a model based on SUSYSO( 10) x f/( 2 ) using a 
126-dimensional Higgs[31]. The CKM matrix of the quark mixings defined 
by Vqkm = V u lVJ l , i s usually parametrised by [15] 


Vckm ~ 


/ 1 - A 2 /2 

-A 

y AA 3 (1 — p — irj) 


A AA 3 (p — ip) 
1 - A 2 /2 AA 2 

-AA 2 1 


(33) 


where A = 0.22 and |A| = 0.90. For our choice of the diagonal up-quark mass 
matrix in Eq.(7), we have V u l = 1 leading to Vckm = Vj L . Since rrid = 
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we have V eL = Vj L 
have 


Vckm■ Neglecting higher power of A in Eq.(33), we 


V f 


eL 



(34) 


which is almost same as V e L given in Eq.(32) except the difference in sign 
before some entries. The positions of the zeros are the same. Such linkage 
gives partial justification to our motivation for the choice of the charged 
lepton mass matrix (15). 


4 Summary and Discussion 

The left-handed Major ana neutrino mass matrix uill which explains the 
inverted hierarchical pattern of neutrino masses, has been generated from the 
seesaw mechanism using non-diagonal texture of the right-handed Majorana 
neutrino mass matrix Mr and diagonal form of the Dirac neutrino mass 
matrix . We have explained the leptonic mixing matrix generated from the 
diagonalisation of rn^L of the inverted hierarachical model and the mixing 
angles so far obtained sm 2 2# 12 ~ 0.999, is too large for the explanation of 
the LMA MSW solution. Such high value of sin 2 26 V2 can be tonned down 
by considering the contribution from the charged lepton mass matrix having 
special structure in the 1-2 block. With such consideration, the predictions on 
the mixing angles at the high scale are sin 2 2# 12 = 0.8517, sm 2 20 23 = 0.9494 
and V e3 = 0.159 which are consistent with the LMA MSW solution. 

The above results which are calculated at the high energy scale (say, 
ju = M u = 2 x 10 16 GeE) where the seesaw mechanism operates, decrease 
with the decrease in energy scale, under the quantum radiative corrections 
within the framework of the MSSM. This is a good feature at least for sin 2 2 #12 
in this inverted hierarchical model as it gives the stability under quantum 
radiative corrections and shows complete consistency of the model with the 
explanation of the LMA MSW solution. The present finding fails to support 
the claim that an arbitrary mixing at the high scale can get “magnified” 
to a large mixing, and even possibly maximal mixing at the low scale [28]. 
Experimental data from a Neutrino factory may confirm the pattern of the 
neutrino masses in near future, and hence the sign of Am| 3 . Such confirma¬ 
tion of the detailed pattern of neutrino masses and their mixing angles is very 
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important as it may give a clue to the understanding of quark masses and 
their mixing angles within the framework of an all-encompassing theory[5]. 

Though we have constructed both uill and m e in a model independent 
way and have shown how the inverted hierarchical model of neutrinos can 
explain the present experimental data particularly LMA MSW solution, the 
present work is expected to be an important clue for building models from 
the grand unified theories with the chiral U(l) symmetry. 
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